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Abstract 

We point out a generic inconsistency of the coupling of ordinary gravity as de- 
scribed by general relativity with matter invariant only under unimodular diffeo- 
morphisms (TDiffs), and some previously studied exceptions are pointed out. The 
most general Lagrangian invariant under TDiff up to dimension five operators is 
determined, and consistency with existing observations is studied in some cases. 
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1 Introduction 



There is a well known way of obtaining the General Relativity Lagrangian which is associ- 
ated w,th the „a,„e of Fey„,„a„ 18|. although ma.y othe. scie.t.t. have eo.t.buted to ,t. 
starting with Kraichnan [TT] q The idea is the following: if one starts from the Fierz-Pauli 
Lagrangian (which describes free spin two particles in Minkowski space) 

where all indices are raised and lowered with the flat Minkowski metric; in particular: 

h = t^Kp (2) 

it so happens that the equations of motion D^^ = \^§^ are transverse, i.e. 

d'-D^J = (3) 

In order to coupling the graviton field hf^,y to an scalar field (p, say, it is natural to try the 
coupling to the conserved energy-momentum tensor (suitably symmetrized if needed, for 
example using the Belinfante technique), that is 

Li = W^'T,, (4) 

But when this term is added to the matter Lagrangian in an freely falling inertial frame 

Ll = \f''d^<\>dA - ^(0) (5) 

the former energy-momentum tensor is no longer conserved, and the gravitational equation 
of motion is inconsistent. This leads to a series of modifications that eventually end up in 
the Hilbert Lagrangian. The quickest path to is probably Deser's, [Tj using a first order 
formalism. The aim of the present paper is to explore what room is left in this argument 



^Somc further references can be found in the review article 2^ or in the book by Ortin 12J 



for less symmetric nonlinear completions, notably the ones we dubbed TDiff, which are 
invariant under coordinate transformations whose jacobian enjoys unit determinant. These 
have been explored in [1], where further references can be found. 



2 The linear approximation 

It is nevertheless clear that given a consistent theory (such as General Relativity itself) its 
linear part in any analytic expansion should be consistent as well (up to linear order). The 
object of our concern in the present paper will be the linear deviations from flat Minkowski 
space, i.e. 

g^lv = Vi^u + K,hf,^ (6) 

where 7/^,^ is the Minkowski metric, and = SttG. This equation is taken to be an exact 
one; it can be looked at as the definition of h^^. 

Now, it is a fact of life that Ip = i^\J^ has got dimensions of length, and that Mp = 
enjoys dimensions of mass. The value of Newton' s constant indicates that at the scale of 
terrestial experiments, Mp ~ lO^^GeV. This means that the field /i^^^ enjoys the proper 
canonical dimension (one) of a four- dimensional gauge field. 

The inverse metric is defined as a formal power series: 

Diffeomorphisms with infinitesimal parameter ^'^ act on the full metric as 

SQi^u = £{09^iu (8) 

whereas in terms of the fluctuations 

SK. = edpKu + - {d^^u + d,^^) + h^^d.C + Kud^C (9) 
This is, again, an exact formula, in the sense that there are no k corrections to it. 



The symmetry as above, without any restrictions, is the one corresponding to General 
Relativity, and in the present paper will be referred to simply as Diff. When the vector 
is restricted to 

a„r = (10) 

the symmetry is broken to what we call T(ransverse)Diff [6]. 
The total action is then defined by 

S=^SH + Sm (11) 

Here Sh represents the purely gravitational sector, which is the most general Lorentz in- 
variant dimension four operator that can be written with the field hais, and its derivatives. 
It can be parametrized by a string of constants, namely the ones associated with the kinetic 
energy, Cj,z = 1 . . .3 and the ones associated to the potential energy for the fiuctuations, 
which is the most general quartic potential in the fiuctuations hap, namely, Aj, i = 1 ... 11. 
The overall scale of the potential energy is related to the cosmological constant, A = M^. 
Cosmological observations seem to favor the tiny value Md ~ 10~^eV. 

M^(l + ^A,.+ ^(A./.^-2A3/.„,^"0 + 

{X,h' + 8A5 V/i'^^/i^ - GXehh^ph'^^) + 

[Mh" - UX^h^K^h'^^ + 32A9/i V/i^^/^^- 
ASX^oKph^^h^sh'^ + 12An [Kph^^f^)) (12) 

Let us remark, first of all, that the structure of the symmetry transformations of both 
Diff and TDiff is such that terms in the Lagrangian of 0(Mp") are related to terns of 
both 0(Mp") as well as, because of the abelian part, terms of 0{Mp^^^). This means 
that in the variations of the kinetic energy part we can keep only the piece in Mp , since 
the other part (0(Mp)) of the variation should cancel with the Mp contribution to the 



kinetic operators of order O(jg^), which we have not considered. The only piece we can 
consistently consider there is then the Fierz-Pauli abelian part given by 

5 V = Mp (a^e. + (13) 

The situation is different, however, in the potential energy piece. In order to cancel the 
Fierz-Pauli variation of the C'(Mp^) term, it is necessary to consider the 0{Mp) variation 
of the 0{j^) term. This means that the full action is invariant under the full variation 
([9]) up to dimension five operators (which means 0{jj^) in the kinetic energy part, and 
Oi^jprs) in the potential energy piece). 
Under those provisos, TDiff needs that 

Ai = A3 = As = Aio 

A2 = Ag = Ag = All 

A4 = As (14) 

The most general TDiff invariant potential depends on four arbitrary parameters. In 
some studies it is frequent to restrict the gravitational equations of motion to the linear 
approximation; this means quadratic terms in the gravitational lagrangian. From a field 
theoretical viewpoint there is no reason lo leave away any relevant (in the renormalization 
group sense) operators. 

• Let us first consider the case Aj = Vz. This corresponds to vanishing cosmological 
constant in general relativity. First of all, TDiff enforces Ci = 1. 

Besides, there are two exceptional values, namely Cj = when TDiff is enhanced 
to full Diff . This is the only combination for which the wave operator is transverse. 

+ jDhr]o.p (15) 
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Indeed 

= ^^Dd^Kf, + ^i^dpdpd^h'"^ + ^^^Dd^h (16) 
By the way, it is worth noticing that the metric condition 



^^9al3 = (17) 

is identically satisfied to o{k) and poses no restriction on /i^,^. 

The other remarkable value is ci = 1 , C2 = | , C3 = | where the symmetry is enhanced 
with a Weyl invariance, denoted by WTDiff, and the wave operator is traceless in 
the absence of cosmological constant. To be specific, 

r/^^^;. = (c3 - ^) Oh + (C2 - I) d.d.h'^'' (18) 

The analysis in [1] shows that these two are the only instances where only spin two 
is present, with no scalar contamination. 

Let us now consider the effect of Aj 7^ 0. First of all, Diff invariance is recovered 
when Aj = 1 Vi. Curiously enough, as such, and in the quadratic approximation, 
the term 

< - -^^^ (19) 

p 



as well as 



- (20) 



do have the interpretation of masses. 

Only when the background around which we perturb is not flat, but a constant curva- 
ture space, with metric ^^jy, do these paramerers recover the meaning of cosmological 
constant. In that case it is mandatory to substitute all derivatives by background 
covariant derivatives, i.e., 

dahf,u -> 'Vahf.u (21) 
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and to raise and lower indices using the corresponding background metric: 

^ g,a-up^^^ (22) 

Let us study now consistency of the couphng, which was the main motivation of 
this work. The matter Lagrangian has been denoted by 5'^- Up to dimension five 
operators for a scalar field, which in a free falling locally inertial reference system has 
lagrangian 

Ll = lv"'d,<Pd.<P - V{<P) (23) 

assuming a Z2 symmetry 

0^-0 (24) 

the allowed matter operators when a gravitational field is present can be parameter- 
ized by three constants, /ii . . . fXs: 

(25) 

Remember that the variation of a scalar field is 

H^i^d,(t> (26) 

In order to enjoy TDiff invariance, it is necessary that yUi = 1. Diff invariance needs 
in addition that ^2 — — ^- The matter equations of motion are 



5Sm 



^ _n0 _ v'{<j>) + ^ (^/.la. {h^^dp<j>) - ^d^ {hd^^) - /i3^V^'(0)) (27) 



5(t) 

and the gravitational equations 

5{^ + Sj)_. 1 /I . ^ , 1 /„ ,,,,, /"2,^,,2^ 



There is generically no problem of consistency, except in the two exceptional cases. 

First of all, when Sh has extended Diff symmetry, and consequently a transverse wave 
operator, this forces to have the same Diff symmetry through the linear expansion 
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of the \/\g\ term; otherwise consistency of the couphng enforces extra condictions 
on the matter, a very weird situation indeed (i.e. Bianchi identities are still valid on 
the gravitational side, so that by consistency the same identities must hold true on 
the matter side as well). Nevertheless, it is not fully devoid of interest to study the 
situations in which there are exceptions to this rule; this we did in a previous work 

When there is WTDiff symmetry, it is clear that the matter lagrangian should also 
be scale invariant in order for the corresponding energy-momentum to be traceless. 
In our example, this corresponds to //i = 2^2 and /X3 = 0. 

There are models in which Diff invariance in the matter sector is reached using in the 
volume element some other scalar density, such as the square root of the determinant 
of a matrix built out of fields and their derivatives (as in the very interesting ones 
proposed in [9]) instead of the ^/\g\ term implicit in the metric volume element. The 
tensor that appears as the source of gravity in Einstein's equations is covariantly 
conserved thanks to the equations of motion of the fields in the scalar densitjo. Of 
course that tensor is not the usual energy-momentum tensor of General Relativity, 
which now is not conserved. The reason is that in order for the Rosenfeld energy- 
momentum tensor to be equivalent to the canonical Belinfante one what is needed is 
not only Diff invariance, but also the standard metric volume element [5j. This topic 
seems worthy of some further investigation. 



3 Observational constraints 

In this section we will outline the way in which one can constraint the space of paramaters 
of the linearized theory (i.e. Cj, and Aj) using experimental results on deviations from 
•^Although its flat limit seems to be different from the canonical energy-momentum tensor. 
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Newton's inverse square law. For simplicity we will illustrate with a very particular example 
so no definite conclusions can be drawn concerning the viability of this kind of models. 

Detailed computations of the propagators can be found in [4j, where the authors con- 
sidered a gravitational Lagrangian f|T2|) with all Aj = except = ^^^2 which has 
the interpretation of a mass for the scalar part of the graviton, present generically in this 
kind of models with TDiff invariance. It turns out that for a conserved energy-momentum 
tensor coupled to gravity in the form 

Li = i V i^iT'"' + ^^2V'"'T) (29) 



then in momentum space the interaction is 
Li = 

where we have defined 



K2 H — 1^1 -ITT^ 2 (30) 



A;2 V 2 / AcP 



1 3 

Ac = C3 - - + C2 - -cl (31) 

with the constraint Ac < because of unitarity [4J. The first term corresponds to the 
usual spin 2 exchange while the second one is an additional massive scalar interaction. 
Let us turn our attention to a particular example, namely the matter Lagrangian (l25il . 
Unfortunately the corresponding energy momentum tensor is not conserved. However, a 
conserved tensor can be defined as 

= r,.-^^..(^^(5p0)' + (/i3-l)m)+O(^) = 

]-d,<t>dA - lv,u (kd,<pr - V{<P)] + 0(-V) (32) 



In the particular case that fi^ = 2^2 ~ 1 (which includes the DifT invariant Lagrangian) our 
energy-momentum tensor can be written in terms of the new one and its trace in such a 
way that the coupling M~^h^i,T^^'^ is of the form ( l29l) with 

2 

(33) 



Now we can apply directly the preceeding results and study experimental constraints to this 
model. The exchange of additional massive scalar degrees of freedom produces a Yukawa 
like potential which is usually parametrized as [Ij 

y(r) ~ ^ (l + ae-^) (34) 

The parameter a is then the ratio between the spin 2 and the scalar couplings, in our 
particular case 



" = ^-72 = T^T— (35) 



Ac kI a Ac 

While A gives the range of the interaction, or equvalently the mass of the scalar exchanged 



= ^ (36) 

Notice that one has to impose m| < since as we have said absence of ghosts requires 
Ac < 0. 

There are important constraints on the strength of hypothetical Yukawa interactions 
for a wide range of A. Through (|35|) and fl36|) it is then possible to constraint the space 
of parameters of the linearized theory. We will use figures 4, 5 and 9 of reference [T], 
which show regions allowed and excluded for a corresponding to A in the ranges 10~^m- 
10~^m, 10~^m-10~^m and 10~^m-10^^m respectively. Since we are just interested in general 
behaviours and not in accurate results we will approximate the experimental curves by 
straight lines. The original plots are in logarithmic scale so we have experimentally allowed 
regions of the form 

|a|</cA" (37) 

We just have to substitute this expresion into to get bounds for our parameters. There 
are however four parameters to play with (/i2, rn"^, C2 and C3). First, it is interesting to 
see the order of magnitude for the mass once we fix the values of C2 and C3. The result 
is ploted in FiglH It can be seen that greater values for the mass are favoured, being the 
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|X2 VS. IrPg^l 




Figure 1: The shadowed region shows experimentally allowed values for |m|| (in m ^) and /U2 
for given values of C2 and C3, expressed in terms of Ac, and in the range A G 10^^m-10~^m. 

lower bound around Irn^] ~ 5 x lO^^m^^ ~ 0.02 eV^, and that the allowed region rapidly 
decreases with |Ac|. 

Another possibility is to fix and fi2 and see, in the plane (02,03), how far from 
Diff invariance (which corresponds to C2 = C3 = 1) can we move away. Remember that 
we also have to take into account the restriction Ac < 0. For the first range, A G 10~^m- 
lO^^m there is no hope of seeing an experimental curve that appreciably deviates from the 
parabola Ac = because of ( l36l) and the tiny values of A. Increasing Img] does increase 
the allowed region, which is between both curves, but does not produce a plot in which 
the curves are visibly separate. It can be understood if we realise that increasing the mass 
also increases C3 on the parabola through C2 and (l35i) - (!36l) . An approximate definition of 
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the separation could be 
Sep ~ 



(38) 



{C'i)par I - C2(A,m2) + |C2(A, 771,2)2 

where {c^)par and {c^)cur means the value of C3 on the parabola and the experimental curve 
respectively. For the first experimental range one has Sep << 1 independently of the mass 
and in the whole interval. The other two cases do not have that property, which is of 
course related also with the particular values of k and a in flH7|) . Examples of resulting 
plots are FigsEllS] and HI 

Once again the experiment prefers greater values for the mass. All the plots have /i2 = 0, 
other values just move the experimental curve along the parabola, but the qualitative result 
remains unchanged. 



C3 vs. C2 
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Figure 2: Experimentally allowed region in the plane (02,03), for a couple of values of the mass, 
in the range A G lO^^m-lO^^m. The plot is restricted to the zone where the curve appreciably 
deviates from the parabola Ac = 0. 
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C3 VS. C2 
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Figure 3: Experimentally allowed region in the plane (02,03), for a couple of values of the 
mass, in the range A G 10~^m-10^^m. We only show the positive C2 branch. The parabola is 
indistinguisable from the C2 axis. 



C3 vs. C2 
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Figure 4: Experimentally allowed region in the plane (02,03) and in the range A G lO^^m-lO^^m 
for a very tiny mass. That allows us to see the parabola, which was hidden in the previous figure. 
We only show the zone closest to the 03 axis. 
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4 Conclusions 



In this paper we have studied at the hnearized level the viability of gravity models with 
a restricted symmetry, both from the theoretical and observational points of view. While 
the existing observational constraints on additional Yukawa like gravitatory interactions 
do not seem to be a major obstacle, a consistency problem has been identified. At the non 
linear level it appeared as an integrability condition on Einstein's equations Here we 
turned our attention to the linear level in order to see if the problem could be avoided, 
and if so in what type of more clever non linear completions. The main conclusion is that 
it is not generically possible to couple matter (i.e., with an arbitrary equation of state) to 
gravitation in such a way that this coupling has a restricted symmetry only (what has been 
called TDiff ) whereas the purely gravitational sector enjoys a higher symmetry, namely the 
standard Diff invariance, or an additional Weyl symmetry (WTDiff). That this is possible 
in some restricted cases has been already found in a previous paper [3]. 

The condition for arbitrary TDiff matter to be able to couple to Diff gravity without 
restrictions can be stated somewhat more formally by saying that the Rosenfeld (metric) 
energy-momentum tensor has got to be equivalent to the Belinfante canonical form. 

On the other hand, there is a widespread urban legend asserting that unimodular theo- 
ries are equivalent to General Relativity with a cosmological constant. Specific calculations 
both here and in our previous paper [3] have proven it to be groundless. It is a fact that 
in some TDiff models there is no exponential expansion at all which is well known to be 
the benchmark of a (positive) cosmological constant in General Relativity. Therefore, that 
models provide a counterexample to the statement above. 

Nevertheless, as with all legends, there is some partial truth in it. The equations of 
motion of the example in [3] correspond to ci = C2 = C3 = //i = 1 and fj.2 = fis = K = 0, 
that is 

= -^d^^d.c^ (39) 
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whereas the hnear equations of General Relativity with cosmological constant read: 

+ ^Va, = ^ (d^<Pd,<f> + \ - \{d,<j>f^ r^^,^ (40) 

Now, the equations (IHUj) are inconsistent as such, in the sense that only a subsector of the 
theory, namely, the one that obeys 

V-\{d,<pf = C (41) 

can be coupled to gravitation. There are many sectors of matter in a freely falling inertial 
system that do not obey thisurestriction. Actually, together with energy conservation, the 
aforementioned equation implies that both the kinetic and potential energy ought to be 
constant: 

21^(0) = E + C 

(a,0)2 = E-C (42) 

It is clear that for a scalar field in flat space most initial conditions lead to configurations 
that violate those equations. This would mean that an inconsistency would show up once 
a gravitational field is turned on, however weak. More formally, something very strange 
should happen when changing the reference frame from an inertial (freely falling one) to 
another in which a gravitational field is present. 

Those are the reasons what we say that the coupling is generically inconsistent. Let 
us accept nevertheless, for the sake of the argument, that physics is so restricted. Then a 
glance at the equations (HOl) of General Relativity shows that they are indeed equivalent 
to flHUj) provided we identify 

But this is only due to our choice of the arbitrary constants, and under the asumption that 
the coupled sector is only the one that obeys (HI]), a deeply misterious condition from a 



*This physically means that the pressure vanishes, (for a perfect fluid the generally covariant lagrangian 
can be identified with the physical pressure cf.jT^), i.e. that the matter is what cosmologists call dust 
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General Relativistic perspective. In the general TDiff case the analogous condition to (HTj) 
is 

- ^^{d,^f + (-"3 - 1)^(0) = C (44) 
and the system is not equivalent to General Relativity with a cosmological constant. 
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